Quantum size effects on spin-tunneling time in a magnetic resonant tunneling diode 
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We study theoretically the quantum size effects of a magnetic resonant tunneling diode (RTD) 
with a (Zn,Mn)Se dilute magnetic semiconductor layer on the spin-tunneling time and the spin 
polarization of the electrons. The results show that the spin-tunneling times may oscillate and a 
great difference between the tunneling time of the electrons with opposite spin directions can be 
obtained depending on the system parameters. We also study the effect of structural asymmetry 
which is related to the difference in the thickness of the nonmagnetic layers. It is found that the 
structural asymmetry can greatly affect the traversal time and the spin polarization of the electrons 
tunneling through the magnetic RTD. The results indicate that, by choosing suitable values for the 
thickness of the layers, one can design a high speed and perfect spin-filter diode. 
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I. INTRODUCTION 

The field of semiconductor spintronics has attracted a 
great deal of attention during the past decade because 
of its potential applications in new generations of tran- 
sistors, lasers, and integrated magnetic sensors. In ad- 
dition, magnetic resonant tunneling diodes (RTDs) can 
also help us to more deeply understand the role of spin 
degree of freedom of the tunneling electron and the quan- 
tum size effects on spin transport processes [l|. By em- 
ploying such a magnetic RTD, an effective injection of 
spin-polarized electrons into nonmagnetic semiconduc- 
tors (NMSs) can be demonstrated I* 1 this regard, 
the II- VI diluted magnetic semiconductors (DMSs) [J, Q 
are known to be good candidates for effective spin injec- 
tion into a NMS because their spin polarization is nearly 
100% and their conductivity is comparable to that of typ- 
ical NMS. A very promising II- VI DMS for spin injection 
is (Zn,Mn)Se, which has been previously used for spin 
injection experiments into GaAs @ and ZnSe 0- The 
(Zn,Mn)Se-based RTD with highly spin-polarized elec- 
tron current has been suggested by Egues § and exper- 
imentally demonstrated by Slobodskyy et al. Also, 
different types of ma gne tic RTDs have been proposed 
both theoreticall vJ]^JOJ]lJp,[ll[l5l[l^| and experi- 
mentally m H MSa El El . 

One of the key parameters in operation of magnetic 
RTDs is the time aspect of tunneling process, which has 
been the focus of much research in the past decade, be- 
cause it is an important parameter for better understand- 
ing of the spin-dependent tunneling phenomena in high- 
speed devices. This quantity may strongly depend on 
the quantum size of the devices, however, to the best of 
our knowledge, no theoretical study on the dependence of 
quantum size on the spin-tunneling time in the magnetic 
RTDs has so far been reported. 
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Recently, based on the group velocity concept, several 
theoretical studies of tunneling time in different mag- 
netic junctions have been done. Guo et al. found ob- 
vious features of separation of spin-tunneling time in 
ZnSe/Zn^Mn^Se heterostructures [HI, [H, [25[. Zhai 
et al. [261 ] studied the tunneling time in magnetic bar- 
rier structures consisting of two identical or unidenti- 
cal magnetic barriers and magnetic wells. Wang et al. 
[27| investigated tunneling properties of spin-polarized 
electrons traversing ferromagnetic/insulator (semicon- 
ductor) double junctions and reported that the tunneling 
time strongly depends on the spin orientation of tun- 
neling electrons. The effects of Rashba spin-orbit in- 
teraction [H, [2t| and Dresselhaus spin-orbit coupling 
|30| on the traversal time of ferromagnetic/ semiconduc- 
tor/ferromagnetic heterostructures have also been inves- 
tigated. By considering the Rashba spin-orbit coupling 
in the semiconductor and significant quantum size simul- 
taneously, it has been found that, as the length of the 
semiconductor increases, the spin-tunneling time will in- 
crease with a behavior of slight oscillation, whether for 
the spin-up electrons or for the spin-down ones [29] . Fur- 
thermore, it has been demonstrated that the Dresselhaus 
spin-orbit coupling, unlike the Rashba spin-orbit interac- 
tion that damps the motion of electrons, does not prolong 
the traversal time of electrons tunneling through the het- 
erostructures [3fj| . 

The aim of the present study is to investigate the 
quantum size effects of a typical magnetic RTD on spin- 
tunneling time and its dependence on the concentration 
of magnetic ions. Our device, similar to that used in 
Ref. [9j , is based on a quantum well made of diluted mag- 
netic semiconductor (Zn,Mn)Se between two (Zn,Be)Se 
barriers and surrounded by highly n-type ZnSe layers. In 
such a structure, the energy levels of the quantum well 
states depend on spin direction due to the exchange split- 
ting, meaning that the energy levels of the (Zn,Mn)Sc 
layer for spin-up electrons will be different from that of 
spin-down electrons. We will show here that this spin 
splitting of the energy levels, which is controlled by an 
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applied magnetic field and the value of Mn concentra- 
tion, enables one to select the resonant condition for the 
desired spin by adjusting the quantum size of the de- 
vice. This paper is organized as follows. In Sec. II, we 
present the model and formalism for traversal time of 
electrons through ZnSe/ZnBeSe/ZnMnSe/ZnBeSe/ZnSe 
magnetic RTD. Numerical results and discussions for 
spin-tunneling time and the degree of electron-spin polar- 
ization in both the symmetric and the asymmetric struc- 
tures are presented in Sec. III. We conclude our findings 
in Sec. IV. 



II. MODEL AND FORMALISM 

Consider a spin unpolarized electron current injected 
into a ZnSe/ZnBeSe/Znx^Mn^Se/ZnBeSe/ZnSe struc- 
ture shown in Fig. 1, in the presence of magnetic and 
electric fields along the growth direction (taken as the 
z-axis). The conduction electrons that contribute to the 
electric currents interact with the 3d 5 electrons of the Mn 
ions via the sp-d exchange interaction. Hence, the exter- 
nal magnetic field, B, gives rise to the spin splitting of the 
conduction band states in the Zni_ x Mn x Se layer. There- 
fore, the injected electrons see a spin-dependent poten- 
tial. Due to the absence of any kind of scattering center 
for the electrons, the motion along the z-axis is decou- 
pled from that of the x — y plane, which is quantized in 
the Landau levels with energies E n = (n+ ^)huj c , where 
n = 0, 1, 2, • • • and hcu c — eB/m* . In such a case, the 
motion of electrons along the z-axis can be reduced to 
the following one-dimensional Schrodinger equation 



h 2 d 2 ^ az {z) 
2m* dz 2 



+ U„ m (z)tli* m (z)=E z tli„ m (z) , (1) 



where the electron effective mass m* is assumed to be 
identical in all the layers, E z is the longitudinal energy 
of electrons, U az (z) is the effective potential seen by a tra- 
verse electron and is given as U Gz (z) = V s + Uq — eV a z/L 
in the ZnBeSe layers (0 < z < L\ and L\ + L 2 < z < 
L) where Uq is the height of the ZnBeSe barriers and 
U<y z {z) = V S + V x {z) + V Ux {z) - eVaz/L in the ZnMnSe 
layer {L\ < z < L\ + L 2 ). Here, L\ and L s are, respec- 
tively, the widths of left and right ZnBeSe layers, and L2 
is the width of the ZnMnSe layer (L = L\ + L 2 + £3); 
V s = \g s P"B G • B describes the Zeeman splitting of the 
conduction electrons, where a is the conventional Pauli 
spin operator; V x (z) is the heterostructure potential or 
the conduction band offset in the absence of a magnetic 
field, which depends on the Mn concentration x and is 
the difference between the conduction band edge of the 
ZnMnSe layer and that of the ZnSe layer; V az (z) is the 
sp-d exchange interaction between the injected electron 
and the Mn ions and can be calculated within the mean 
field approximation. Hence, the sum of the last two terms 
can be written as 

V x {z) + V as (z) 



lAE(x) N aa zXcS SB s (j^^ 
xe(z-L 1 )e(L 1 +L 2 -z) , 



where 



AE(x) = E g (x) - 
= -0.63X 



E g (0) 

\- 22x 2 - 195x 3 + 645a: 4 



(2) 



(3) 



is the sum of the conduction and valance band offset 
under zero magnetic field, when the real (effective) Mn 
concentration is x (x e ff = x[l — x} 12 ). Here, Bs(- ■ •) is 
the Brillouin function and S = | is the spin of the Mn 
ions. a z — ±i (or f , J.) are the electron-spin components 
along the magnetic field. We should note that this form 
of Eq. [3]is valid only for < x < 0.1 31]. The last term 
in U az (z) denotes the effect of an applied bias V a along 
the z axis on the system. 

In order to study the tunneling time of electrons 
through the structure, we adopt the group velocity ap- 
proach [32I [33] ]. in which the tunneling time of a spin- 
polarized electron can be defined as r CT= = j dz/v gMz (z), 
where the spin-dependent group velocity, v g ^ z , is defined 
as the ratio of the average probability current density 
S az = Ke\h(ij) a -^dit'a-, /dz)/im*] to the probability den- 
sity li/vj 2 of the particle |23l. [24. \2l\ . In this regard, the 
spin-tunneling time can be written as 



in 



1 



^ Jo l7<r* (2)|Im[tan0<^ ( z )] 



dz 



where 



tan# a;i (z) 



1 dipu^ (z) 



ip (Tz (z) dz 



1<Jz {z) = l -^/2m*[E z -U^(z)} 



(4) 



(5) 



(6) 



In the above equations, ^ az (z) is the spin-dependent 
wave function of the heterostructure. Since we have con- 
sidered that the electrons tunnel through the magnetic 
structure from the left (z < 0) to the right (z > L), 
under the influence of the applied voltage V a , the wave 
functions in each region can be written as 



VV 3 {z) = 

e ik laz z _|_ r ^^ e -ik laz z^ z<0, 

A 2az Ai[p az (z)} + B 2 aBi[p az (z)}, < z < ii, 

A 3<7z Ai[p az (z)} + B 3a Bi[p az (z)} 1 Li <z< Li + L 2 , 

A 4az Ai[p az (z)] + B 4 a z Bi[p az (z)} 1 Li + L 2 < z < L, 



to 



z > L. 



(7) 



Here, 



and t„ are the reflection and the trans- 



mission amplitudes; k\„ z = yj2m*(E z — V s )/Ti and 
k 5lJz — y / 2m*(E z —V s + eV a )/h are the electron mo- 
menta; Ai[p (Ti (z)] and Bi[p (Ti (z)] are the Airy func- 
tions with p (j (z) = [E z — U a (z)]L/(eV a X) and A = 
[-Ti 2 L/(2m*eV a )} 1 / 3 . 
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FIG. 1: (Color online) Spin-dependent effective potential pro- 
file in the magnetic RTD under influence of both the magnetic 
and electric fields along the z axis. The up and down arrows 
represent spin-up and spin-down subbands, respectively. The 
two ZnBeSe layers have the same widths L\ = L3=50 A. 
L 2 =100 A, a;=0.04, B=2 T and V a =5 mV. 



The constants Aj at and Bj az (with j'=2-4) can be de- 
termined from a system of equations formed by ij) <Tz and 
its derivative for the same z value. Then, we can relate 
ip az and its derivative ip' a at two positions z and z' by 
the transfer matrix (M) as follows: 



</v, 0) 



Mn M12 
M 2i M 22 



(8) 



Applying the results of Eq. [8] in Eq. [5l the value of 
tan^ (z), which determines the tunneling time at posi- 
tion z, can be written in terms of its value at position z' 
as 



tan^ (z) 



1 MnlT^^OItan^^O+M 



12 



(z) | Mai ha. (*') |tan0 CT2 (*') + M 



22 



(9) 

Therefore, using Eqs. [9]and|4j the spin-tunneling time 
can be calculated for the desired magnetic RTD. 



III. RESULTS AND DISCUSSION 

In this section, we use the formulas given above 
to investigate the quantum size effect on the spin- 
tunneling time and electron-spin polarization in the 
ZnSe/ZnBeSe/Zni-^Mn^Se/ZnBeSe/ZnSe heterostruc- 
tures. In the numerical calculations we have taken the 
following values: U = 92 meV T = 2.2 K, T = 1.4 
K, and N a = -0.27 eV [Hj], Ep = 5 meV, B = 2 T, 
V a = 5 mV, <7 S = 1.1, and m* = 0.16 m e (m e is the 
mass of the free electron). In the tunneling process at 
low temperatures, the electrons with energy near Fermi 
energy (Ep) carry most of the current; for this reason 
we have done our numerical calculations at E z = Ep. 
Also, due to the band-gap bowing of the Zni-^Mn^Se 
layer [3l|, [34[ , we have examined the effects of three val- 
ues of Mn concentrations (x — 0.02, 0.04, and 0.07) for 
the paramagnetic layer. 
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FIG. 2: (Color online) Spin-tunneling time for electrons 
traversing the symmetric structure: L\ = Lz = 50 A. The 
dashed and solid lines correspond to spin-up and spin-down 
electrons, respectively. 



First, we study the tunneling time in the symmetric 
and the asymmetric structures, depending on the thick- 
ness of the ZnBeSe layers. The structure is called sym- 
metric (asymmetric) if L\ = L3 [L\ ^ L3). Figure 2 
shows the spin-tunneling time as a function of thickness 
L2 of the Zni-zMn-jSe layer. When L 2 is zero, the tun- 
neling time is independent of the spin orientation and 
hence r-f = rj.. At x = 0.02 [Fig. 2(a)] and with increas- 
ing L 2 , the tunneling time oscillates for both spin-up and 
spin-down electrons; however, the length period of oscil- 
lations for Tj is shorter than Tp At this concentration, 
the paramagnetic layer has its minimum value of band- 
gap, which is smaller than the band gap of the ZnSe layers 
[3l| . In the case of x = 0.04 [Fig. 2(b)], the oscillation 
does not change considerably for spin-down electrons, but 
the length period significantly increases for the spin-up 
ones. At x = 0.07, the tunneling time oscillates only for 
spin-down electrons, as shown in Fig. 2(c). It is impor- 
tant to note that at x ~ 0.04 and under zero magnetic 
field, the band gap of Zni-^Mn^Se is nearly the same as 
that of ZnSe, while for larger values of x such as x = 0.07, 
the Zni-^Mn^Se layer behaves as a potential barrier in 
comparison with the ZnSe layers [34| . 

The appearance of oscillatory and non-oscillatory be- 
haviors in the tunneling time is due to the effective po- 
tential U az (2) and reflects this fact that the difference 
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FIG. 3: (Color online) Spin-tunneling time for electrons 
traversing the asymmetric structure: L\ = 10 A and L-j, = 90 
A. The dashed and solid lines correspond to spin-up and spin- 
down electrons, respectively. 
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FIG. 4: (Color online) Spin-tunneling time for electrons 
traversing the asymmetric structure: L\ = 90 A and L3 = 10 
A. The dashed and solid lines correspond to spin-up and spin- 
down electrons, respectively. 



between the position of the Fermi energy and the bottom 
of the conduction band on the energy axis for each spin 
orientation in the paramagnetic layer has dominant effect 
on the behavior of spin-tunneling time. The physical ori- 
gin of the oscillatory behavior is explained by the quan- 
tum well states. As is well known for resonant tunneling 
through double-barrier structures [H, [36| , when the in- 
cident energy of electrons coincides with the energy of a 
quasibound state in the quantum well, a resonance condi- 
tion is fulfilled and the transmission coefficient of the elec- 
trons through the heterostructure strongly increases. On 
the other hand, the position of the quantum well states, 
formed in the Zni-^Mn^Se layer, strongly depends on 
the well thickness £2- Therefore, with continuous varia- 
tion in I/2, the position of the resonant states varies and 
this leads to the oscillations of the tunneling time. As a 
remarkable feature in the oscillations, one can see sharp 
dips in comparison with the broad peaks in the tunnel- 
ing time curves, which correspond to the narrowing of 
the width of resonant bands arising from confinement of 
electrons to the paramagnetic layer. With increasing L 2 , 
these narrow bands (levels) quickly cross the energy of in- 
cident electrons in the left ZnSe layer on the energy axis, 
the tunneling time decreases, and the sharp dips appear. 
In contrast, the peaks are broad, due to the broad gaps 



between the discrete resonant levels in the ZnMnSe layer. 

Figs. 3 and 4 show the spin-tunneling times in two 
asymmetric structures: (i) case L\ < L3 where L\ = 10 
A and L3 = 90 A; (ii) case L\ > L3 where L\ = 90 
A and L3 = 10 A. In the case of L\ < L3 (Fig. 3), 
the spin-tunneling times are almost increasing functions 
of L2 and there is no oscillation between r-f and rj, in 
all three Mn concentrations. In contrast, for the case of 
L\ > L3 (Fig. 4), the spin-tunneling times show oscil- 
latory behavior very similar to those of the symmetric 
structure [L\ = L3). The reason for discrepancy in the 
tunneling time of two asymmetric structures is that the 
voltage drop within the ZnBeSe and ZnMnSe layers de- 
pends on the position z [see the last term in U az (z)], and 
hence, the depth of quantum well in the paramagnetic 
layer strongly depends on the position of the layer with 
respect to the origin. Therefore, in the case of L\ > L3, 
the incident electrons see a deeper quantum well with 
respect to the the case of L\ < L3, and this affects the 
features of the resonant states. If we increase the bias 
voltage, the oscillatory behavior in the asymmetric struc- 
tures with Li < L 3 appears, too. These results may be 
important from experimental point of view and indicate 
that special care must be taken during sample growth in 
order to make a magnetic RTD with low power consump- 
tion, high speed, and greater spin-filter efficiency [bj ] . It 
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with increasing L2, the spin polarization retains positive 
and its value increases almost linearly and reaches nearly 
100%, indicating an excellent spin filtering effect. This 
means that for this value of x and in all the structures, 
the tunneling process of spin-up (spin-down) electrons is 
always a slow (quick) process. In the cases of x = 0.02 
and x = 0.04, however, the spin polarization can change 
sign for particular ranges of Li in the structures with 
Li = L 3 and L\ > L 3 . Therefore, the case of r-f > 
or Tj < Tj may strongly depend on the thickness of the 
paramagnetic layer. 

According to the above results, the tunneling process of 
the spin-polarized electrons through the magnetic RTD 
can be divided into slow and quick processes. However, 
we cannot say which one of the spin orientations of the 
tunneling electrons always corresponds to the slow pro- 
cess and which one corresponds to the quick process. 
Such a feature occurs for Mn concentration x < 0.04 
in both the symmetric and the asymmetric (L\ > L3) 
structures. We would like to point out here that our 
calculations have been performed under the assumption 
of a phase-coherent tunneling process, which applies to 
heterostructures with narrow wells and barriers. When 
the heterostructures become thicker, we should replace 
the phase-coherent tunneling by a sequential process and 
these features will change. 



FIG. 5: (Color online) Degree of spin polarization as a func- 
tion of the Zni-^Mn^Se layer thickness L2 for (a) symmetric 
and [(b) and (c)] asymmetric structures. 



IV. CONCLUSION 



is necessary to point out that the effects of the width 
of the paramagnetic layer on the spin-dependent current 
densities in such a magnetic RTD have been studied and 
the oscillatory behaviors in the current-voltage charac- 
teristics reported [Til ]. 

Now, for further understanding of the quantum size 
effects of the system on the spin-tunneling times, we cal- 
culate the degree of spin polarization of the tunneling 
electrons, which can be defined as P = (r-f — rj.)/ (r-f +tj_). 
The results are plotted in Fig. 5 for both the symmetric 
and the asymmetric structures. It is clear that the value 
of spin polarization can be greatly changed by the Mn 
concentration, the thickness of the ZnMnSe layer as well 
as the status of structural symmetry. At x = 0.07 and 



Using the group velocity concept and the particle cur- 
rent conservation principle, we have shown how the ge- 
ometry and the size of the device affect the spin-tunncling 
time and the degree of spin polarization of tunneling elec- 
trons in a (Zn,Mn)Se-based magnetic RTD. The tunnel- 
ing time for spin-up and spin-down electrons and hence 
the degree of spin polarization may strongly depend on 
the width of the paramagnetic layer and the Mn concen- 
tration. We found that, due to the oscillatory behavior 
of spin-tunneling time with increasing thickness of the 
(Zn,Mn)Se layer, special care should be taken for design- 
ing a magnetic RTD with high efficiency Furthermore, 
the present results may open a new way to control the 
degree of spin polarization and design the high speed 
magnetic devices. 
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